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ABSTRACT 

The focus of this work was to explore the numerical simulation of the 
probabilistic aerothermodynamic response of propulsion system components to 
randomness in their environment. The reusable rocket engine turbopumps were 
selected as an example because of the severe cryogenic environment in which 
they operate. The thermal and combustion instabilities, coupled with the 
engine thrust requirements from start up to shut down, lead to randomness in 
the flow variables and uncertainties in the aerodynamic loading. 

The probabilistic modeling of the turbopumps aerodynamic response was 
accomplished using the panel method coupled with Fast Probability 
Integration methods. The aerodynamic response in the form of probabilistic 
rotor blades and splitter loading were predicted and the results presented 
for specified flow coefficient and rotor preswirl variance. Possible future 
applications of the aerothermodynamic probabilistic modeling in engine 
transient simulation, condition monitoring and engine life prediction are 
briefly discussed. 


*Work funded by Space Act Agreement C-99066-G. 



INTRODUCTION 


Traditional engineering design and/or analysis procedures are 
deterministic in nature. The uncertainties in the load definition, material 
properties and component characteristics are handled in these designs 
through safety factors. Recently, a major effort was initiated by the NASA 
Lewis Research Center under the PSAM (Probabilistic Structural Analysis 
Methods) program to develop analytical and numerical techniques for 
predicting the probabilistic response of the structural components [1). 

This new approach makes it now possible to assess the risk of failure 
inherent In a given design and to identify the design features affecting 
risk [2]. Examples of the application of this new methodology to propulsion 
system components can be found in reference [3], 

The purpose of this work is to explore the use of probabilistic 
modeling techniques in predicting the aerothermodynamic response of 
propulsion system components to randomness in their environment. The 
randomness In the flow variables and the associated uncertainties in a 
component's aerothermodynamic performance can be a consequence of combustion 
instabilities, thermal cycling and/or the engine built variance. A 
probabilistic flow analysis was used to predict the aerodynamic performance 
parameters In the form of distribution functions, for specified variance in 
the component's environment. The analysis was also used to Investigate the 
sensitivity of the aerodynamic response to the various parameters. The 
results obtained from the probabilistic determination of the SSME high 
pressure turbopump performance are presented, and their possible use in 
condition monitoring, design optimization and ultimately life prediction are 
discussed. 
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OVERVIEW OF PROBABILISTIC SIMULATION METHODOLOGY 


The purpose of probabilistic modeling is to simulate the effect of 
variability In a propulsion system environment on the performance of its 
components. The random variables are those with a certain amount of 
statistical uncertainty. In the case of turbomachines, they might Include 
the incoming flow conditions, the rotational speed, the turbomachinery blade 
configuration or other parameters which can affect the performance. The 
simulation of the corresponding variation in the performance requires an 
analysis based on probabilities. The aerodynamic response of the system can 
be expressed as a function of the random variables as follows: 

Y - f(x) (!) 

where Y Is the vector of response variables and x Is the vector of random 
variables. Random variables may be characterized experimentally or 
theoretically in terms of distribution models. Experimental measurements 
from a number of tests may give a set of data points which are scattered 
over a certain range. The characterization of these variables based on 
experimental results can be depicted in a histogram, in which each block 
size represents the number of occurrences within a small range. The 
probability density function, or the smooth curve representing the scatter 
in the data, can be obtained from a sufficiently large number of tests. 

This can easily be simulated using a statistical package, to determine the 
best analytical fit for data distribution (normal, log normal, Weibul, 
etc.). For example the best statistical representation for the data of 
Fig. 1 representing the measured SSME high pressure turbopump RPM at 100% 
RPL (rated power level) , was found to be a normal distribution with a mean 
value of 34,133 and a standard deviation of 210. 
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The response variable Y, which is a function of random variables, will 
have its own probability density and cumulative distribution functions. 
Numerical methods are often used to obtain solutions to equation (1) , since 
the multiple integrals in the solution can be evaluated analytically only if 
the function f(x) is linear, the random variables are normal and their 
number is restricted. In addition, the functional relationship of equation 
(1) is usually not given in a closed form but rather at a discrete number of 
points from a deterministic flow solution. 

Several computational approaches have been considered to evaluate the 
probabilistic response in the PSAM program including the Monte Carlo 
Simulation, the method of moments and the fast probability integration 
methods. A discussion of the capabilities of these methods and examples 
involving comparison of the results of their applications to sample problems 
with a small number of random variables, can be found in, reference [4]. The 
FPI fast probability integration code developed under PSAM was used in the 
present study, since it is generally applicable when the relationship 
between the response variable and the primary variables is not given in a 
closed form. It is based on constructing an analytical approximation 
(polynomial) for the performance function and using a modified Newton 
iteration scheme to compute the perturbed solutions about a mean state [4] . 
This approach was demonstrated to have high accuracy with respect to the 
Monte Carlo simulation and required less computational effort [4] , and to 
exhibit good convergence properties for small changes about the 
deterministic solution at the expansion point [5]. 
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PROCEDURE 


For a given propulsion system or component, it is necessary to first 
identify the uncertainties in its environment that can affect its 
aerothermodynamic performance. The probabilistic aerodynamic response to 
the associated variance in the appropriate parameters is then computed. The 
determination of the probabilistic aerodynamic response requires 
probabilistic tools as well as a deterministic aerodynamics code. The 
appropriate choice of the latter will depend on the particular component, 
the simulated flow phenomena and the aerodynamic response parameters of 
interest. It can be a viscous or inviscid code for unsteady or steady, 
three dimensional or other reduced levels of approximating the flow field in 
the propulsion system. The FPI method was used in the probabilistic 
analysis of the deterministic solutions obtained for the different specified 
values of the random variables. 

SSME TURBOPUMFS 

The SSME reusable high performance rocket engines operate in a severe 
environment over a range of variable thrust. The engine thrust requirements 
from start-up to full power steady loading is followed by throttling to SO- 
TO percent thrust and then to full power before shut down, all occurring 
within about 450 seconds (Fig. 2). The thermal loading and combustion 
instabilities in the complex flow environment lead to randomness in the 
performance parameters and uncertainties in the component's loading. Since 
the turbopumps of the SSME operate in environments that are more severe and 
uncertain than in other noncryogenic turbomachinery applications, they were 
identified as primary candidates for developing probabilistic aerodynamic 
models. The probabilistic results are presented for the high pressure fuel 


5 



turbopump. The same methodology can be used to predict the probabilistic 
aerodynamic response of other propulsion system components after the proper 
identification of the sources of randomness in their environment and of the 
appropriate analytical and numerical methods for each particular 
application. 

OPERATING ENVIRONMENT 

The SSHE is controlled by an on-board computer which monitors the 
engine's vital functions, performs engine throttling and initiates shutdown 
procedures in the event that redline parameters are violated. The engine 
data handled are identified by Pi's (parameter identification number) and 
stored in the Command and Data Simulator System (CADS) . Figures 3a and 3b 
from reference [6] show the variation In some of the engine's performance 
parameters with the power level over the 0.65 to 1.09 RPL range. The mean 
value of the various parameters can be seen clearly in the figures, as well 
as an Indication of their variance. The variance in some of the performance 
parameters at the 0.65 and 1.09 RPL Is better illustrated in Figs. 4 through 
9. These test data can be used to generate probability density functions 
and to determine the type of variance and standard deviation for the 
corresponding parameters using statistical tools. The mean values of the 
flow inlet parameters and the pump's rotational speed were taken according 
to Rocketdyne's power balance model [7]. 

The random variables, those with uncertainties that can affect the 
turbopump performance include the flow inlet conditions (P ,T t ), the mass 

flow rate and the rotational speed. For simplicity, the geometric 
configuration was not Included in the present study. When expressed in 
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dimensionless parameters, the performance of turbomachines can generally be 
shown to be dependent on a reduced number of variables. In the case of 
incompressible flow in the pump rotor, two parameters affect the performance 
of its rotor, namely the flow coefficient, C^/U (which is directly related 

to the mass flow), and the flow preswirl parameter, C^/U. A simple 

derivation based on the turbine's Euler equation can show the linear 
relationship between the specific work or head and these two parameters. 

The response variables presented in the results are the nondimens ional axial 
and circumferential blade forces. 

RESULTS 

Figure 10 shows the high pressure turbopump blade geometry at the tip, 
with three splitter blades in each blade- to-blade passage. The presented 
results for the flow field were obtained using the PANEL dode [8] because of 
its ability to handle multibody configurations with good resolution at the 
blade leading edges. 

Table 1 lists the performance parameters at 90% RPL according to 
reference [7], These were used as the mean values in the present 
probabilistic flow analysis. The computed blade surface pressure 
distribution at the mean values of the parameters are presented in Figs. 11 
through 14 . 

The probabilistic flow analysis results were obtained using the FPI 
code. They are presented in the form of cumulative distribution functions, 
CDF, for the dimensionless axial and tangential blade loading, in response 
to specified variance in the flow coefficient, C^/U, and the preswirl 

parameter, C /U. The sensitivity, S, of the computed blade loading to the 
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two primary random variables (C^/U, C^/U) were also computed. They are 

identified as S(C^/V) and S(CyU) in Figs. 15 through 24, which also shows 

the cumulative distribution function, CDF, for the blade loading. Two types 
of flow coefficient and preswirl variance were considered a normal 
distribution and a log norm (Table 2) . Figures 15 through 18 present the 
results obtained for the normal distribution with 10% standard deviation. 

The computed probabilistic aerodynamic results for a log normal C^/U, C u /U 

distribution with 5% standard deviation are presented in Figs. 19 through 
22. Some of the results for a normal distribution with 5% standard 
deviation are presented in Figs. 23 and 24. Comparing Figs. 19a with 23 and 
22 with 24, one can see that both the aerodynamic response and its 
sensitivity to C^/U, C^/U are almost the same for the log norm and normal 

distribution. Table 3 summarizes the probabilistic aerodynamic response 
results of Figs. 15 through 24. The sensitivity analysis results shown in 
the figures indicate that the blade loading sensitivities to the preswirl 
C^/U and to the flow coefficient C^/U are of the same order of magnitude. 

As the blade loading increases, the sensitivity to C^/U increases and to 

C x /U decreases in most cases. The probabilistic axial loading of the first 

blade and its sensitivities need further investigation because of the 
difficulties encountered in its FPI analysis. 

DISCUSSIONS 

The effect of changing the flow coefficient and preswirl on the blades 
aerodynamic loadings are listed in Table 4 and presented in Figs. 25 and 26. 
These figures confirm the linear relationship between the blade loading and 
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both random variables over a range of ±10% variation in C x /U and C^/U . This 

makes it possible to evaluate analytically the probability distribution 
function of the response variables for the case of normal distribution of 
the variance in the two random variables. Such a procedure is recommended 
for assessing the accuracy and evaluating the range of perturbations over 
which the FPI results are within the desired accuracy. In general, more 
than two random (primitive) variables affect the steady state performance of 
turbomachines in the case of compressible flow. The engine transient 
performance is affected by an even larger number of random variables [9] in 
which case the analysis can become more complex, and the probabilistic 
results more difficult to assess. 

CONCLUSIONS 

One possible application of the presented turbopump's probabilistic 
aerodynamic response is to assess the uncertainties in the SSME Internal 
axial loads. According to reference [10], the large axial loads (100's k 
lbf levels) mandates axial load analyses at all SSME operating modes to 
ensure the adequacy of the balancing piston and bearing designs. This will 
require an integration of the axial blade loadings over each blade surface 
area at the different operating conditions. While the presented 
dimensionless blade loading results only depend on two parameters, the 
actual blade loading is affected by the dynamic head (based on the relative 
velocity) and the blade area. Additional uncertainties can be introduced by 
these parameters as can be seen from the measured variance In the inlet 
total flow properties of Figs. 5, 6, 8 and 9. Such a probabilistic axial 
loading analysis can be accomplished by expanding the scope of the present 
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exploratory investigation with more computational efforts to simulate the 
probabilistic blade loading distribution over the SSME's operating range. 

RECOMMENDATIONS 

In general, the results of a probabilistic aerodynamic analysis can be 
used in the definition of the aerodynamic loads in structural analysis. For 
such an application, It is recommended that additional parameters, such as 
variance in the individual blade geometries, be included In the simulation. 
Presently, the probabilistic simulation already accomplished in the PSAM 
program include material properties, blade characteristics and mechanical 
loading, but the aerodynamic loads have not been defined probabilistically. 
However, under steady state conditions, the aerodynamic loadings are 
generally very small compared to centrifugal loading and to the thermal 
stresses in Internally cooled blades [11). The author's own experience [9 ] 
suggests that the probabilistic modeling of the system's transient 
aerodynamic response would produce more dramatic results. Such an analysis 
can be very useful for predicting probabilistic engine responses including 
turbopumps, surge and cavitation based on measured or deduced variance in 
the important performance parameters. The results of a probabilistic 
aero thermodynamic analysis can also be very useful in failure and life 
prediction [12]. In particular, the thermal shock transients in the rocket 
engines can have a strong Impact on the life of its components. Results 
showing the complex steady state aero thermodynamic flow fields in the SSME 
turbines were presented In reference [13]. A probabilistic 

aero thermodynamic analysis would provide valuable information on the thermal 
blade loading and its impact on life. Condition monitoring is an area in 
which the probabilistic aerothermodynamic modeling can find applications. A 
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conveniently measureable flow property can be identified as a monitoring 
candidate based on the sensitivity results of the probabilistic aerodynamics 
analysis. A probabilistic analysis will make it possible to take into 
consideration the variance in the monitored response with engine built. 
Finally, design optimization studies can be based on the probabilistic 
aerodynamic results, with a better definition of the uncertainties than in 
current deterministic methods. 
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SSME TURBOPUMP'S PERFORMANCE DATA AT 90% RPL 
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TABLE 2 

RANDOM VARIABLES 
HPFP AT 90% RPL 


C /U 
x' 



Distribution 

Mean 

Value 

Standard 

Deviation 

Distribution 

Mean 

Value 

Standard 

Deviation 

1 

normal 

0.2467 

10% 

normal 

0.293 

10% 

2 

log normal 

0.2467 

5% 

log normal 

0.293 

5% 

3 

normal 

0.2467 

5% 

normal 

0.293 

5% 
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TABLE 3 


PROBABILISTIC AERODYNAMIC RESPONSE OF 
THE SSME HPTP BLADES AT THE TIP 


Blade 

Loading 

Case 

Mean 

Value 

1 

Standard 

Deviation 

Case 

Mean 

Value 

2 

Standard 

Deviation 

Case 

Mean 

Value 

3 

Standard 

Deviation 

F ti 

0.5292 

0.0306 

0.5275 

0.01595 

0.5275 

0.01595 

F al 

-- 


0.8945 

0.0971 

0.8945 

0.0971 

F t2 

0.4610 

0.0506 

0.4565 

0.0263 

0.4565 

0.0263 

F a2 

0.5175 

0.0550 

0.5666 

0.0292 

0.5666 

0.0292 

F t3 

0.4320 

0.0454 

0.4275 

0.02375 

0.4275 

0.02375 

F a3 

0.6545 

0.0541 

0.6500 

0.0283 

0.6500 

0.0283 

F t4 

0.4985 

0.0584 

0.4918 

0.0303 

0.4918 

0.0303 

F a4 

0.63175 

0.0663 

0.6255 

0.0346 

0.6255 

0.0346 


Cases correspond to those of Table 2. 




TABLE 4 

EFFECT OF C /U AND C /U ON BLADE LOADING 
x u' 



C /U 
x' 

C /U 
u 

F tl 

F al 

F t2 

F a2 

F t3 

F a3 

F t4 

F a4 

1 

0.2467 

0.293 

0.528 

0.898 

0.455 

0.566 

0.427 

0.649 

0.491 

0.624 

2 

0.222 

0.293 

0.501 

0.889 

0.413 

0.521 

0.39 

0.607 

0.441 

0.568 

3 

0.2344 

0.293 

0.517 

0.898 

0.435 

0.544 

0.409 

0.63 

0.4666 

0.597 

4 

0.259 

0.293 

0,538 

0.894 

0.476 

0.588 

0.445 

0.668 

0.516 

0.652 

5 

0.2714 

0.293 

0.554 

0.9005 

0.503 

0.617 

0.4685 

0.6952 

0.5475 

0.6805 

6 

0.2467 

0.264 

0.541 

0.936 

0.458 

0.569 

0.427 

0.654 

0.492 

0.6275 

7 

0.2467 

0.278 

0.516 

0.887 

0.4395 

0.5474 

0.412 

0.629 

0.474 

0.604 

8 

0.2467 

0.308 

0.54 

0.906 

0.4725 

0.5858 

0.443 

0.671 

0.509 

0.646 

9 

0.2467 

0.322 

0.514 

0.854 

0.456 

0.563 

0.428 

0.645 

0.493 

0.621 
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90203^ 902038 902042 902046X10 1 

TEST NUMBER 

T I CURE 1. - TEST STAND DATA SHOWING THE SKATTER IN THE 
MEASURED HPFP SPEED WITH ENGINE BUILD (RET. 14). 
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POWER LEVEL IN PERCENT. RPL 



TIME. SEC 


FIGURE 2. - TYPICAL SPACE SHUTTLE MAIN ENGINE 104 PERCENT 
MISSION. 


274X10 2 



902037 902039 902041 902043 902045 902047 

TEST NUMBER 

FIGURE 4. - MEASURED HPFP SPEED AT 65 PERCENT RPL. 


375x1 0 2 



(a) HPFT SPEED VARIATION WITH POWER LEVEL (REF. 6). 


75x I0 2 



POWER LEVEL 


(b) HPFP DISCHARGE PRESSURE VARIATION WITH POWER LEVEL 
(REF. 6). 

FIGURE 3. - ENGINE PERFORMANCE VARIATION. SAMPLE STANDARD 
DEVIATION. 2.00 SIGMA FROM THE SAMPLE AVERAGE. 



902037 902039 902041 902043 902045 902046X10 1 

TEST NUMBER 

FIGURE 5. - MEASURED HPFP INLET PRESSURE AT 65 PERCENT 
RPL. 
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HPFP IN PRESSURE, TOTAL 



TEST NUMBER 


FIGURE 6. - MEASURED HPFP INLET TEMPERATURE AT 65 PER- 
CENT RPL. 



TEST NUMBER 


FIGURE 8, - MEASURED HPFP INLET PRESSURE AT 109 PER- 
CENT RPL. 


368x1 0 2 



TEST NUMBER 


FIGURE 7. - MEASURED HPFL SPEED AT 109 PERCENT RPL. 



TEST NUMBER 


FIGURE 9. - MEASURED HPFP INLET TEMPERATURE AT 109 PER- 
CENT RPL. 
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PRESSURE COEFFICIENT 









SENSITIVITIES AND CUMULATIVE DISTRIBUTION FUNCTION 



TANGENTIAL FORCE COEFFICIENT, F a2 


(a) TAN6ENTIAL FORCE COEFFICIENT; MEAN VALUE, 0.4610; 
STANDARD DEVIATION, 0.0506. 



AXIAL FORCE COEFFICIENT, F t 2 


(b) AXIAL FORCE COEFFICIENT; MEAN VALUE, 0.5175; STANDARD 
DEVIATION, 0.0550. 

FIGURE 16. - PROBABILISTIC AERODYNAMIC LOADING AND SENSITIV- 
ITIES SECOND BLADE TANGENTIAL FORCE NORMAL C x /U, Cy/U DIS- 
TRIBUTIONS WITH 10 PERCENT STANDARD DEVIATION. 



(a) TANGENTIAL FORCE COEFFICIENT; MEAN VALUE, 0.4320; 
STANDARD DEVIATION, 0.0454. 



.4 .5 .6 .7 .8 .9 1.0 


AXIAL FORCE COEFFICIENT, F a3 

<b) AXIAL FORCE COEFFICIENT; MEAN VALUE, 0.6546; STANDARD 
DEVIATION, 0.0541. 

FIGURE 17. - PROBABILISTIC AERODYNAMIC LOADING AND SEN- 
SITIVITIES THIRD BLADE TANGENTIAL FORCE NORMAL C/U, 
C u /\] DISTRIBUTIONS WITH 10 PERCENT STANDARD DEVIATION. 
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SENSITIVITIES AND CUMULATIVE DISTRIBUTION FUNCTION 



(a) TANGENTIAL FORCE COEFFICIENT; MEAN VALUE, 0.4985; 
STANDARD DEVIATION, 0.0584. 



(b) AXIAL FORCE COEFFICIENT; MEAN VALUE, 0.63175; STANDARD 
DEVIATION, 0.0663. 


FIGURE 18. - PROBABILISTIC AERODYNAMIC LOADING AND SEN- 
SITIVITIES FOURTH BLADE AXIAL FORCE NORMAL C/U, C/U 
DISTRIBUTIONS WITH 10 PERCENT STANDARD DEVIATION. 



TANGENTIAL FORCE COEFFICIENT, F n 


- (a) TANGENTIAL FORCE COEFFICIENT; MEAN VALUE, 0.5275; STAND- 

3 ARD DEVIATION, 0.01595. 

1 





AXIAL FORCE COEFFICIENT, F al 

(b) AXIAL FORCE COEFFICIENT; MEAN VALUE, 0.8945; STANDARD 
DEVIATION, 0.0971. 

FIGURE 19. - PROBABILISTIC AERODYNAMIC LOADING AND SEN- 
SITIVITIES FIRST BLADE AXIAL FORCE LOG NORMAL C x /U, 
Cy/U DISTRIBUTIONS WITH 5 PERCENT STANDARD DEVIATION. 
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SENSITIVITIES AND CUMULATIVE DISTRIBUTION FUNCTION 



(a) TANGENTIAL FORCE COEFFICIENT; MEAN VALUE, 0.4565; 
STANDARD DEVIATION, 0.0263. 




AXIAL FORCE COEFFICIENT, F a2 

<b) AXIAL FORCE COEFFICIENT; MEAN VALUE, 0.5666; STANDARD 
DEVIATION, 0.0292. 


Q 


> 


to 

w 



AXIAL FORCE COEFFICIENT, F a3 

(b) AXIAL FORCE COEFFICIENT; MEAN VALUE, 0.65; STANDARD 
DEVIATION, 0.0283. 


FIGURE 20. - PROBABILISTIC AERODYNAMIC LOADING AND SENSI- 
TIVITIES SECOND BLADE AXIAL FORCE LOG NORMAL C x /U, Cy/U 
DISTRIBUTIONS WITH 5 PERCENT STANDARD DEVIATION. 


FIGURE 21. - PROBABILI I STIC AERODYNAMIC LOADING AND SEN- 
SITIVITIES THIRD BLADE AXIAL FORCE LOG NORMAL C x /U, 
C y /U DISTRIBUTIONS WITH 5 PERCENT STANDARD DEVIATION. 
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SENSITIVITIES AND CUMULATIVE DISTRIBUTION FUNCTION 



(a) TANGENTIAL FORCE COEFFICIENT; MEAN VALUE, 0.4918; 


STANDARD DEVIATION, 0.0303. 



(b) AXIAL FORCE COEFFICIENT; PCAN VALUE, 0.6255; STANDARD 
DEVIATION, 0.0346. 


FIGURE 22. - PROBABILISTIC AERODYNAMIC LOADING AND SENSI- 
TIVITIES FOURTH BLADE AXIAL FORCE LOG NORMAL C/U, C y /U 
DISTRIBUTIONS WITH 5 PERCENT STANDARD DEVIATION. 



FIGURE 23. - PROBABILISTIC AERODYNAMIC LOADING AND SENSI- 
TIVITIES FIRST BLADE TANGENTIAL FORCE NORMAL C/U, C/U 
DISTRIBUTIONS WITH 5 PERCENT STANDARD DEVIATION. 
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AXIAL FORCE COEFFICIENT, F aq 

(b) AXIAL FORCE COEFFICIENT; MEAN VALUE, 0.6255; STANDARD 
DEVIATION, 0.03%. 

FIGURE 24. - PROBABILISTIC AERODYNAMIC LOADING AND SENSI- 
TIVITIES FOURTH BLADE AXIAL FORCE NORMAL C x /U, C^/IS DIS- 
TRIBUTIONS WITH 5 PERCENT STANDARD DEVIATION. 



c x /u 


FIGURE 25. - VARIATION OF BLADE LOADING WITH FLOW CO 
EFFICIENT. 
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